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ABSTRACT 

The Central Molecular Zone (CMZ) of the Milky Way shows several peculiar 
properties: a large star formation rate, some of the most massive young star clusters 
and molecular clouds in the Galaxy, and a twisted ring morphology in molecular 
gas. In this paper, I use SPH simulations to show that most of these properties can 
be explained as due to a recent outburst of AGN activity in Sgr A*, the central 
super massive black hole of the Milky Way. In particular, the narrow ring of dense 
gas, massive gas clouds, young star clusters and an elevated SFR can all be caused 
by the passage of an AGN outflow through the system, which compresses the gas and 
triggers fragmentation. Furthermore, I show that the asymmetric distribution of gas, 
as observed in the CMZ, can be produced by outflow-induced instabilities from an 
initially axisymmetric gas disc. Angular momentum mixing in the disc produces some 
low angular momentum material, which can subsequently feed Sgr A*. These processes 
can occur in any galaxy that experiences an AGN episode, leading to bursts of nuclear 
star formation much stronger than pure bar-driven mass inflows would predict. 

Key words: quasars: general — accretion, accretion discs — ISM: evolution — 
Galaxy: centre — stars: formation — open clusters and associations: general 


1 INTRODUCTION 

Over the past two decades, various observations revealed 
that most galaxies harbour super massive blac k holes 
(SMBHs) at their ce ntres ([Kormendv &; Richstonel Il995l : 
iMagorrian et al.l fl998). It is expected that these SMBHs 
and their immediate surroundings spend a few percent of 
their lifetimes as active galactic nuclei (AGNs), accreting 
gas at rates exceeding a few percent of their Eddington 
limit MEdd = 2.2 x 10 -8 (Mbh/M©)M© yr -1 . During this 
time, the AGN luminosity can exceed that of the whole 
host galaxy. The AGN radiation can heat the surround¬ 
ing gas by photoionization and Comp ton heating. More 
impo rtantly, quasi-relativistic winds (Tombes^e^ffi 010b, 
20141) and massive mole cular outflows ([Sturm et al 2011 : 
Cicone et al.l 120141 . l2015h have been detected in AGN, th e 
latter most likely driven by the former (Tombesi et al. 2015). 
The observed outflow energy rate is ~ 5% of the AGN lu - 
minosity, as pre dicted by analytical (IZubovas &; King|l2012h 
and numerical ([Costa et al.l 12014 120151 ) models. Such out¬ 
flows affect the AGN host galaxy on scales ranging from the 
immediate AGN environment to the circumgalactic medium. 
Importantly, footprints of these effects can remain visible 
for a much longer time than the duration of the AGN phase 


itself, allowing us to probe the past activity of presently 
quiescent galactic nuclei. 

A prototypical and well-studied quiescent galactic nu¬ 
cleus is our own Galctic Centre (GC), which contains a ra¬ 
dio sour ce Sgr A*, identi fi ed as an Mbh = 4.3 x 10 6 M© 
SMBH ([Ghez et al.l 120081 : iGillessen et al.l 120091 ) . The cur¬ 
rent bolometric luminosity of Sgr A* is v ery low, Lboi ~ 
300L© ~ 10 -9 LEdd (e.g., iMelia &; Falckel l200lh . There is 
evidence that this low luminosity is caused by two effects: 
there is little gas supplie d into the several-p c wide sphere of 
influe nce around Sgr A* (|Bower et al.l [20031 : iMarrone et al.l 
2006), and the gas that is su pplied is accreted with very low 
efficiency ([Wang et al .1120131 ). 

However, there is growing evidence that Sgr A* used 
to be much brighter in the past. Light echoes from molec¬ 
ular clouds in the Galactic centre reveal that Sgr A* X- 
ray l uminosity was Lx ^ 10 39 erg s" 1 about 100 years 
ago dRevnivtsev et al.l |2004 iTerrier et al.l l2010h and per- 
haps on several o t her o ccasions in the past ~ 10 3 years 
(|Gando Rvu et all l2012h . On longer timescales, there is 
evidence of significant activity some time during the 
past 10 Myr. This evidence comes from the presence of 
young stars in the centr al parsec ([Paumard et al.l [2006; 
iHobbs &; Navakshinl l2009h and the two large gamma-ray 
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emitting bubbles (the Fermi bubbles: ISu et al.l 2010b which 
have properties consistent with being inflated se veral Myr 
ago by either a n AGN iet (I Guo &; Mat hewsll20 12h or wind- 
driven outflow ([Zubovas et al.1 12011 : IZubovas &; Navakshinl 
12011 also see Section [2] below). In particular, in these two 
papers my colleagues and I showed that an Eddington- 
limited AGN outburst that started 6 Myr ago and lasted for 
1 Myr would have produced two bubbles with sizes, shapes 
and energy content consistent with observations. 

Given this evidence of increased past activity, it is inter¬ 
esting to check whether footprints of that activity might be 
present closer to the GC, such as in the Central Molecular 
Zone (CMZ) of the Galaxy. The C MZ is a ring-like structure 
of predominantly molecular gas (IMorris &; Serabvnl Il996l : 

I Jones et"ahl l201ll : iMolinari et al.l l201ll ). It is oriented par¬ 
allel to the plane of the Galaxy, with a radius Rcmz ~ 
200 p c and thickness hcmz ~ 100 pc dMorris &; Serabvnl 
I1996T ). giving an approximate aspect ratio H/R ~ 0.25. 
The total gas mass contained in this region is Mg mz ~ 
10 8 Mfi (IMorris &; Serabvnl Il996l : IPierce-Price et all l2000l : 
IMolinari et al.ll2011 ). 

The CMZ contains several star formation regions, such 
as the molecular cloud Sgr B2 (the mo st massive GMC 
in th e Galaxy; Ms gr B 2 — 2 x 10 6 M iRe vnivtse v et al.l 
120041 ) . Sgr C and Sgr D ([Pierce-Price et ah 2000h . Two 
young massive s tar clusters, Arches ( t age — 3 d= 1 Myr, M g 
4 — 7 x 10 4 M(^. iMartins et al.l 120081 : IPortegies Zwart et al.l 
120021 : iFiger et al.l 2002| ) and Q uintuplet (t age ~ 4 ± 1 Myr, 
M cx 1O 4 M 0 , Figer et ahlll999 ) are also associated with the 


region. The mean sta r formation rate in the region, M* ^ 
0.06 — 0.08 M 0 yr -1 (lYusef-Zadeh et al.ll2009l : llmmer et al.l 
; alF 


120121 : iLongmore et al.l 120131 ). would deplete the molecu¬ 


lar gas reservoir in tde 


10 yr < tGaiaxy ~ 10 yr, 


therefore it seems likely that the present star formation 
rate is elevated by some recent process, i.e. one that oc- 
curred _at most a few dy namica l times ago (see, however, 
iRodriguez-Fernandez fe Combesl 120081 . who show that the 
elevated star formation rate can be maintained by material 
transport via bars from further out in the MW). 

Recent Herschel observations reveal a complex struc¬ 
ture of num erous clo uds and filaments comprising the CMZ 
(IMolinari et al.ll201ll ). They are mostly positioned on a nar¬ 
row elliptical ring, which is twisted and offset from the dy¬ 
namical centre of the Galaxy. The line-of-sight velocities of 
the clouds are larger than would be expected if the gas was 
moving along the ring. In particular, the “20 km/s” and “50 
km/s” clouds certainly do not follow the ring. Therefore the 
ring is most likely just a transient structure, created by an 
accidental accumulation of gas. 

In this paper, I present numerical simulations designed 
to show what effect a recent burst of AGN activity in Sgr A* 
might have had upon the CMZ properties. Starting from ide¬ 
alised initial conditions, the AGN wind shapes the CMZ into 
its present-day elliptical configuration and triggers gravita¬ 
tional instabilities which lead to formation of massive young 
star clusters. The overall conclusion is that AGN activity 
might be as important a factor in CMZ evolution as the 
perturbation due to the bar potential of the Galaxy. 

The paper is structured as follows. I briefly review the 
AGN wind-driven outflow model and derive the important 
quantities regarding its effect upon the cold CMZ gas in 
section [2j Then, I describe the numerical simulation setup 


(section (3} and present their results (section 0|) , paying par¬ 
ticular attention to the morphology of the gas and orbits of 
star clusters, the star formation rate and the global proper¬ 
ties of the star cluster population. I discuss the main impli¬ 
cations of these results in section [5] and summarize them in 
section [6j 


2 AGN WINDS AND THEIR EFFECT ON 
DENSE GAS 


2.1 Wind feedback model 


Observations show that a large fraction of AGN have fast 
(Vw ^ 0.1c} winds with kinetic power Lki n ~ 0.05 Lagn 
(|Tombesi et al.l lOlOal fbh . The ubiquity of observed winds 
suggests that they are emitted with a wide ope ning an¬ 
gle, perhaps quasi-spherically (|Nardini et al.ll2Q15l ). These 
winds shock a gains t the galactic ISM and drive large- 
scale outflows ([Kind 12010 ) . In principle, if the shock oc¬ 
curs very close to the SMBH, the inverse-Compton cool¬ 
ing allows the wind to rapidly lose mos t of its ener gy and 
the resultant flow is momentum-driven (iKin 3(2003). How- 
ever, in practice, the electron-ion thermalisation timescale 
in the two-temperature plasma of the shocked wind is 
long enough so that c ooling is e fficient only at dist ances 
< 1 pc from the AGN (iFaucher-Giguere fc Quataertl 12012 ) 
and the outflows are essentially always energy-driven. 
Nevertheless, in a non-uniform ISM, most of the out¬ 
flow energy leaks out through low-density channels, ef¬ 
fectively re covering the mom e ntum-driven solution for the 
dense gas iNavakshin J2014I : IZubovas &; Navakshinl 120141 : 
iGabor fe Bournaud 20141 ). Such outflows are then good can¬ 
didates for communicating the AGN luminosity to the host 
galax y and establishing the M — a relation ([Kin 3 '2003 . 
l201ft) . The ram pressure of the wind can only push away 
dense gas if the AGN luminosity is higher than the critical 
value L a ^ 3 x 10 45 crf O o er g s _1 , where a m lOOcrioo km/s 
is the velocity dispersion in the host galaxy. Assuming the 
AGN to be Eddington-limited, this translates into a mini¬ 
mum mass M a ~ 2.3 x 10 7 crf O o M 0 . 

In our Galactic centre, one then expects that any out¬ 
flow would affect the CMZ predominantly via the ram pres¬ 
sure of the wind. This pressure exerts a force 


_ L_ < 47tGMsmbh 


- 1.8 x 10 34 dyn, 


(i) 


where I used the value of Sgr A* mass M — 4.3 x 10 6 M 0 
(|Ghez et al.ll2008l : iGillessen et al.ll2009l ). In comparison, the 
weight of the CMZ is 


Wcmz — 


GMcmzM (< Rcmz) Mcmzct 2 


~ 3 x 


Rqmz Rcmz 

10 34 M 8 crioo^oo dyn, 


( 2 ) 


We see that even an Eddington-limited outburst from Sgr A* 
is unable to lift the CMZ and remove it from the Galaxy. 
Note, however, that the difference between the two quanti¬ 
ties is only a factor ~ 2, so even a modest retention of the 
outflow energy (rather than just its momentum) can dis¬ 
perse the CMZ. This may lead to the CMZ being eroded 
over a longer timescale, as its surface material absorbs some 
of the outflow energy, heats up and evaporates. 
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Sgr A* shapes the Central Molecular Zone 3 


2.2 Mixing of CMZ material 

Even though the CMZ is not dispersed, the AGN wind can 
significantly perturb it. The ram pressure is not negligible 
compared to the weight of the whole CMZ, so a significant 
fraction of the CMZ material is pushed outwards and can 
mix with the rest. 

The CMZ is located outside the sphere of influence of 
the SMBH, so the circular velocity of its components is ap¬ 
proximately independent of radius. Hence the specific an¬ 
gular momentum of CMZ gas increases linearly with radius 
and spans more than an order of magnitude. When gas with 
different angular momentum is mixed together, the result is 
an averaging of this quantity (jHobbs et al.l [2011 b Gas in¬ 
jected from the centre gains angular momentum, while gas 
from outside loses it, in a process similar to a time-inverted 
viscous ring evolution (except that the process responsible 
for the mixing of angular momentum is dynamical perturba¬ 
tions rather than viscosity). We can make a rough estimate 
of the radius at which the CMZ gas would collect after this 
perturbation. If we assume the CMZ follows an RT 2 radial 
density profile with constant H/R , the mass at each radius 
is just M(R —► R + A R) = Mcmz AR/Rcmz- If the CMZ is 
affected out to radius R ou t < Rcmz and all gas is moving 
on circular orbits initially, the total angular momentum in 
the affected region is M ou t'CcircRout/2, therefore the gas cir¬ 
cularizes at a mean radius of Rout/2- If the whole CMZ was 
affected, the circularization radius is 100 pc. Furthermore, 
gas that was initially outside and inside the circularization 
radius moves inward and outward, respectively, increasing 
the density contrast between the gas collected at the circu¬ 
larization radius and gas outside of it. This effect can turn an 
initially radially widespread gas distribution into a narrow 
ring. 

The background density outside the sphere of influence 
of Sgr A*, assuming a singular isothermal sphere potential, 
is 

a 2 

PSIS = 2-xGR 1 - 2,5 X 10_21cr ioo-Rio 2 o g cm -3 , (3) 

which, assuming fully ionised gas of Solar metallicity (/i = 
0.63), gives 

n S is - 2.5 x 10 3 a?oo^roo cm -3 . (4) 

The mean gas density within the CMZ, again assuming an 
R~ 2 density profile, is pcmz ~ 10 -21 MsRiqq g cm -3 ~ 
0.4psis* When the perturbed gas begins circularizing at 
around Rcmz/2, it can easily pass the psis threshold and 
become self-gravitating. 

Gravitational instabilities can lead to formation of fil¬ 
aments and clumps of gas, as well as star clusters. This is 
a stochastic process, so I turn to numerical simulations to 
investigate the properties of the structures forming in a per¬ 
turbed CMZ. 


3 NUMERICAL MODEL 

I run SPH simulations of the CMZ affected by a burst 
of AGN activity. The code used is GADGET-3 , an up¬ 
dated version of the publicly available GADGET-2 ([Springeil 
120051 5. It is a hybrid N-body/SPH code with individual par¬ 
ticle timesteps and adaptive smoothing. 


With these simulations, I intend to show how AGN ac¬ 
tivity can create an asymmetric distribution of gas and star 
clusters from initially smooth and regular conditions. There¬ 
fore, I set up the CMZ as a smooth disc with total mass 
M — 10 8 M 0 , extending from an inner radius n n = 5 pc 
to outer radius r ou t = 200 pc, with a radial density pro¬ 
file p oc r~ 2 . The CMZ has a constant rotational veloc¬ 
ity Vtft = 141 km/s, which makes it rotationally supported 
against the gravity of the background potential. The disc gas 
is initialized with a constant temperature T — 3.2 x 10 4 K, 
giving a constant H/R — 0.25, i.e. the scale height of the 
disc is 50 pc from the midplane at the outer edge. In the ver¬ 
tical direction, the gas density decreases exponentially away 
from the midplane. This set up results in a surface density 
S ^ 1.5 x 10 4 (r/n n ) _1 M 0 pc -2 and a radius-independent 
Toomre Q parameter Q cs 2.5. Thus the disc is initially 
marginally stable against axisymmetric perturbations and 
self-gravity. 

The CMZ has N cmz = 5 X 10 5 particles, giving a par¬ 
ticle mass mspH = 200 M 0 and mass resolution mn Yes = 
40ra sp h = 8000 M 0 . This is fine enough to resolve large 
molecular clouds and clusters. The spatial resolution is vary¬ 
ing due to the nature of SPH. In the initial disc, it ranges 
from ~ 1 pc at the inner edge to ~ 13 pc at the outer edge; 
in the densest regions, the spatial resolution goes down to 
the minimum allowed smoothing length of 0.01 pc. 

The disc is surrounded with a halo of diffuse (p g = 
10 -3 psis) gas extending out to 5 kpc. Although a very 
crude approximation, this density approximately corre¬ 
sponds to the observed gas density in the Milky Way 
(Kal berla &; Pe des 2008). The halo gas temperature is set 
at Th = 2.4 x 10 5 K, equal to the virial temperature of the 
surrounding static, spherically symmetric isothermal back¬ 
ground potential with velocity dispersion a = 100 km/s. I 
comment on the effects of the assumption of spherical sym¬ 
metry in Section l5Tl 

A SMBH with mass Mbh = 4.3 x 10 6 M 0 is embed¬ 
ded in the centre of the gas distribution. The sphere of in¬ 
fluence of the SMBH extends only to Rind ~ 1 pc, there¬ 
fore the SMBH gravity does not affect the dynamics of the 
gas. At the start of each simulation, the SMBH is turned 
on as an AGN of luminosity L = IL^dd and radiates at 
this luminosity for a time t q , after which the luminosity 
drops to zero. The AGN effect on the gas is modelled in 
two ways. First of all, the simulation employs a heating¬ 
cooling function appropriate for optically thin gas of So¬ 
lar metallicity exposed to quasar radiation (|Sazonov et al.l 
2005), with a temperature floor of Tfioor = 10 4 K. Secondly, 
AGN wind feedback is mode lled with the help of virtual par¬ 
ticles (|Navakshin et al.1120091 "). These particles are emitted by 
the SMBH every timestep. They have only trace mass, but 
carry momentum p v irt = kmspucr each, where k = 0.1 is 
an arbitrary constant, and energy E v i r t = Pvirtc/2; p v irt also 
sets the number of virtual particles emitted each timestep. 
They move at a constant velocity v w — 0.1 c radially out¬ 
ward from the source. Whenever an SPH particle contains 
a virtual particle within its smoothing kernel, interaction 
between the two begins. The virtual particle gives up a frac¬ 
tion of its momentum and energy to the surrounding SPH 
particles, weighted by the SPH particles’ smoothing kernel. 
The fraction given up in each interaction is set so that it 
takes ~ 10 steps for each virtual particle to lose 99% of 
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Model ID 

l 

*q (Myr) 

Gink (Myr) 

©ink) (Mg yr- 1 ) 

M sink , tot (1O 6 M 0 ) 

a/b 

Control 

- 

0 

- 

0 

0 

- 

Basel 

1 

1 

4.5 

15.0 

26.7 

1.39 

Base2 

1 

1 

4.6 

11.3 

27.6 

1.34 

Base3 

1 

1 

- 

4.9 

29.7 

1.53 

Base4 

1 

1 

3.4 

12.7 

37.8 

1.39 

Base5 

1 

1 

4.0 

14.2 

37.8 

1.43 

Shortl 

1 

0.7 

2.6 

11.8 

28.1 

1.43 

Short2 

1 

0.7 

3.0 

13.6 

28.9 

1.27 

Short3 

1 

0.7 

3.5 

13.2 

24.3 

1.55 

Short4 

1 

0.7 

3.0 

13.2 

27.9 

1.38 

L0.3 

0.3 

1 


1.7 

10.2 

1.16 

TO.3 

1 

0.3 

1 . 1 

3.7 

23.0 

1.09 

T1.5 

1 

1.5 

- 

2.1 

12.6 

6.78* 

L0.667T1.5 

0.667 

1.5 

- 

2.2 

12.9 

4.23* 


Table 1. Parameters of the numerical models and most important results. The first column shows the model ID. The next two columns 
give the parameters: the AGN Eddington ratio, l = T/I/EdcU an d the duration of AGN phase ir Myr, £ q . The final five columns are the 
primary results: time of onset of significant fragmentation (see text), mean fragmentation rate after the onset, total sink particle mass 
at t = 6 Myr, CMZ ring major-to-minor axis ratio at t = 6 Myr. Models Basel to Base5 and models Shortl to Short5 differ in initial 
microscopic particle distribution, but have identical parameters. Numbers labelled with * use a lower density threshold to identify the 
ring due to lack of significant structures. 



-200 -100 0 100 200 

x [pc] 



-3 - 2-10 1 

log £ [g cm' 2 ] 


Figure 1 . Projected gas density of the Control simulation at 
t = 6 Myr. Several spiral patterns form in the disc, but there is 
no collapse or fragmentation. 


its momentum/energy, at which point it is removed from 
the simulation. This method of transferring feedback self- 
consistently accounts for different optical depth of gas in 
different directions and enables the formation of anisotropic 
feedback-inflated bubbles. In all simulations presented here, 


the wind shock is assumed to be adiabatic, therefore all of 
the wind kinetic energy L w = 0.05Lagn is transferred to 
the gas. 

In order to speed up the simulations and for easier track¬ 
ing of the star formation process, we convert gas particles 
into sink particles according to a Jeans’ condition. Whenever 
a gas particle density increases above the critical value 

( 7T&B T \ 3 —2 K 1 0-12^3 -3 / c \ 

pj = U^gJ m - h " L 5 x i° T4gcm ’ (5) 

nj ~ 1.4 x 10 12 Tl cm -3 , ( 6 ) 

where T 4 = T/10 4 K, it is converted into a sink particle of 
the same mass, which subsequently interacts with the other 
particles only via gravity and can accrete other gas or sink 
particles once they come within 10 -4 pc. The critical den¬ 
sity is such that the corresponding Jeans’ mass is equal to 
raspH- In all simulations, most gas particles stay well below 
the density threshold, with density increasing significantly 
only in self-gravitating clumps on sub-parsec scales; reduc¬ 
ing the threshold by several orders of magnitude does not 
significantly affect the results. 

There are fourteen simulations in total; their parame¬ 
ters - l and £ q - and main results are given in Table [I] I run 
one Control simulation with no AGN activity, in order to 
check whether the initial gas distribution is stableto its own 
self-gravity without external perturbations. Five simulations 
(Basel to Base5) have identical parameters, but use stochas¬ 
tically different initial conditions, as the particles comprising 
the CMZ were cut from different parts of a relaxed glass dis¬ 
tribution. The choi ce of l = 1 and t a = 1 Myr i s motivated 
by previous work ([Zubovas fe Navakshinl Eoid h where we 
found that these parameters produce hot gas cavities very 
similar to the observed Fermi bubbles. Four more simula¬ 
tions (Shortl to Short4) also differ only in microscopic par¬ 
ticle distribution and have t q = 0.7 Myr; this choice is mo¬ 
tivated by the results of the Base simulations, as described 
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Figure 2. Projected gas density of fiducial run Basel at early times. Left: cut side view at t — 0.5 Myr, showing the radial structure 
of the CMZ disc. The central regions are pushed together and begin collapsing vertically, exposing parts of the disc further out to AGN 
feedback. Middle: cut side view at t = t q = 1 Myr. The whole disc is affected, with dense and thick inner ring and surface layers pushing 
into the outer parts. Right: top view at t = t q = 1 Myr. A transition region is visible, with spiral perturbations and self-gravitating 
clumps. 


below. The final four simulations show the effects of varying 
l and tq further, as these parameters are least constrained by 
observations and current models. All but one of these four 
simulations use the same initial conditions as the Basel and 
Shortl models; one simulation, L0.3, crashes before finish¬ 
ing due to one sink particle growing to a mass greater than 
that of Sgr A*, so I reran it with initial conditions identical 
to those of model Base2. All simulations ran for 6 Myr, the 
presumed time interval since the start of the AGN phase. 


4 RESULTS 

I first present the results of the Control simulation, which 
has no AGN activity and is designed to show the stabil¬ 
ity of the initial gas configuration. Then I present the five 
Base models, which all have the same AGN parameters: an 
Eddington-limited outburst lasting for t q = 1 Myr. I discuss 
the morphology and the structures forming in the perturbed 
CMZ, derive fragmentation rates and properties of the ma¬ 
jor sink particle clusters. Later, I comment on the effect of 
varying the AGN outburst duration and luminosity. 

4.1 Control simulation 

In this simulation, the gas disc is allowed to evolve freely 
without external perturbation. Soon after the start of the 
simulation, a circular overdensity forms in the centre of 
the ring and moves outward with velocity v ~ 40 km/s; 
it reaches the outer edge of the disc at t ~ 4 Myr, at which 
point the disc radius has shrunk to ~ 160 pc. This happens 
because the initial conditions feature sharp density cutoffs 
at both the inner and outer edge. As the disc relaxes, density 
waves form, but the one on the outer disc is much weaker 
due to the lower density contrast. 

Subsequently, starting from ~ 1 Myr after the start of 
the simulation, spiral density perturbations form in the inner 



Radius [pc] 

Figure 3. Radial density profiles of CMZ gas in the model Basel 
at t = 0 (black), t = 0.5 Myr (red) and t = 1 Myr (blue). Grey 
points are a sample of SPH particle densities at t = 1 Myr. 


part of the disc. By t = 6 Myr, they are present through¬ 
out the disc, with a stronger density contrast in the inner 
~ 80 pc, reflecting the shorter dynamical timescale there. 
The formation of spiral density waves is expected, since the 
Toomre Q factor of the initial configuration is 2.5, around 
the limit where it becomes unstable to such perturbations 

(Q < 2 ). 

Figure [T| shows the projected gas density map of the 
disc at t = 6 Myr. Despite the presence of both circular and 
spiral density waves, no fragmentation has occurred in the 
CMZ, and no sink particles formed. This lack of fragmen¬ 
tation happens partially due to the imposed temperature 
floor of 10 4 K; if the gas could cool down to much lower 
temperatures, some fragmentation would be expected. How¬ 
ever, stellar feedback and external radiation might keep the 
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Figure 4. Projected gas density of fiducial run Basel after the AGN switches off. All images show top view. Left: t = 3 Myr; spiral 
perturbations are falling inward and facilitating radial migration of CMZ gas. Middle: t = 6 Myr; a narrow asymmetric ring of gas is 
winding around the Galactic centre. White dots show the sink particles; several clusters can be identified. Right: A zoom-in on the most 
massive sink particle cluster. 
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Figure 5. Radial density profiles of CMZ gas in the model Basel 
at t = 1 Myr (black), t = 3 Myr (red) and t = 6 Myr (blue). Grey 
points are a sample of SPH particle densities at t = 6 Myr. 


gas in this configuration above the temperature required for 
fragmentation. I discuss the importance of the gas equation 
of state in Section 15.5.11 


4.2 Fiducial runs: l =* 1, £ q = 1 Myr 

4-2.1 Feedback-driven expansion 

From the start of the simulation, AGN wind impacts the 
inner edge of the CMZ disc and pushes it outward. The ex¬ 
pansion velocity is initially approximately constant, ^ r ,out — 
150 km/s = \/2a. Formally, the purely energy-driven ve¬ 
locity of expansion should be ~ 330 km/s, while purely 
momentum-driven velocity should be ~ 25 km/s. The actual 
velocity is in between the two, suggesting that a large frac¬ 
tion (~ 80%) of the energy is reflected away from the dense 
gas and works to form bubbles in the diffuse halo (this effect 


was investigated in more detail by IZubovas &; Navakshinl 
12014) . Gas temperature at the edges of the CMZ increases 
to > 10 5 K, resulting in some gas evaporating from the sur¬ 
faces and escaping into the diffuse halo; meanwhile, most of 
the gas stays close to the temperature floor. After ~ 0.5 Myr 
(Figure [2] left panel) gas in the inner parts of the CMZ is 
compacted so much that vertical gravitational collapse be¬ 
gins, as the Toomre Q parameter in that region drops to 
< 1. The inner parts of the disc become thinner and less 
exposed to the AGN wind (red line in Figure [3}; meanwhile, 
parts further away become more exposed and start moving 
radially outward. Globally, this results in a transition region 
forming between the evacuated inner part of the disc and the 
stable outer parts. A surface layer of higher density forms 
in the outer part of the disc, where the disc gas is pushed 
inward by the expanding hot bubble in the halo (Figure [2j 
middle panel). 

By the end of the AGN outburst, i.e. at t = 1 Myr, the 
transition region extends from ~ 130 to ~ 160 pc (Figure 
[2] right panel; also blue line and grey points in Figure [3j) . 
Column densities in the transition region exceed 1 g cm -2 ; 
this is enough to bring the Toomre Q parameter below 1 
(Ecrit — OAR^oq g cm -2 ~ 1700R^ 0 qMo pc -2 ). As a result, 
spiral perturbations and self-gravitating clumps form within 
the transition region. 

All five Base models evolve very similarly during this 
stage, since their evolution is dominated by a global process 
(AGN feedback), rather than local gravity. 


4-2.2 Fragmentation 

Once AGN activity ends, the CMZ material collapses back 
toward the centre. The circularization radius is ~ 60 pc, but 
falling material overshoots this limit and moves to radial 
distances of ~ 40 pc before expanding again. In effect, indi¬ 
vidual gas particles move on eccentric orbits and the whole 
disc appears to breathe. The particles on the outer edge of 
the transition region at t — 1 Myr continue moving outward 
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Direction [deg] 

Figure 6. Average radius (black solid line), half-diameter (red 
dashed line) and displacement (blue dot-dashed line) of the dense 
gas ring, defined as n > 100 cm -3 , as function of direction coun¬ 
terclockwise from left in model Basel at £ = 6 Myr. Gas further 
than 200 pc from the centre was excluded when calculating these 
values. The radius (i.e. distance from the origin to the dense ring) 
varies between ~ 80 and ~ 150 pc. The half-diameter (half the 
distance between the ring on opposite sides of the origin) varies 
between ~ 100 and ~ 140 pc, and the SMBH is displaced from 
the centre of the ring by about 20 pc. 

for a while, encompassing the whole disc and causing it to 
expand to a maximum radius of ~ 270 pc at £ — 2 Myr. 

During this time, spiral perturbations continue to grow 
in the disc. By t — 2 Myr, they are clearly visible. By 
t — 3 Myr, they have become significant density enhance¬ 
ments, facilitating infall of material to small radii (Figure 0] 
left panel). The formation and growth of perturbations is a 
stochastic process, and as a result, differences among the five 
Base models begin to manifest. In the rest of this subsec¬ 
tion, the descriptions refer to model Basel, unless otherwise 
noted. 

By the end of the simulation, the CMZ gas is distributed 
strongly anisotropically, with a narrow (d r i ng ~ 20 — 40 pc) 
ring winding around the Galactic centre at distances ranging 
from r m in ~ 80 pc to ~ 150 pc (top-left and bottom in the 
middle panel of Figure [4j respectively; see also Figures [5] and 
ED. The ring has a sharp outer edge, with density dropping 
by several orders of magnitude outside ~ 150 pc, and a less 
pronounced inner edge where density decreases by approx¬ 
imately a factor 4. The ring can be roughly approximated 
as an ellipse with semimajor axis of 140 pc and semiminor 
axis of 100 pc, with the semimajor axis making a 30° angle 
with the Y axis. The SMBH is displaced by ~ 20 pc from 
the geometric centre of the ellipse. The total mass of the gas 
ring, defined as having gas with density ridense > 100 cm -3 , 
is Mdense = 6.9 x 10 7 M©. This value does not include the 
mass of sink particles (see below). The cutoff density is cho¬ 
sen in such a way that if the selected gas cools down isobari- 
cally to 20 K, the column density through the ring would be 
Ydense — Odense X (Xfloor/20K) ^ X Coring — 4 X 10 Cm , 
equal to the threshold used when identifyi ng the ring from 
infrared observations (Mol inari et alJlioilh . 

The densest parts of the perturbations become self- 
gravitating clumps and create sink particles. The first sink 
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Figure 7. Side view of model Basel at £ = 6 Myr. The gas ring 
is slightly disturbed, but generally stays on the Z = 0 plane. 


particles appear before the AGN switches off, but fragmen¬ 
tation starts in earnest at £ > 4.5 Myr. There are five regions 
of intense fragmentation, lighting up one after another. In 
addition, some sink particles form in the dense filaments; 
these most likely represent small star clusters and associ¬ 
ations with masses at or below the simulation resolution 
of 8000M©. By the end of the simulation, the total mass 
in sink particles is M s i n k = 2.8 x 10 7 M©. Thus the total 
mass in sink particles and dense gas is ~ 9.7 x 10 7 M g , i.e. 
almost all of the initial CMZ material ends up either con¬ 
verted into sink particles or in the dense gas ring. The pres¬ 
ence of dense gas is not unexpected, since the initial condi¬ 
tions contain gas denser than the threshold used to identify 
the ring; however, the transformation from a smooth disc 
into a narrow ring is remarkable. The average fragmenta¬ 
tion rate over the 1.5 Myr period of intense fragmentation 
is Mf rag ~ 15 M© yr -1 , but it can rise up to ~ 25 M© yr -1 
(see Figure El- Assuming a typical star formation efficiency 
on such mass scales of 2 — 10% gives a star formation rate 
M* ~ 0.25 — 1.25 M© yr -1 , dropping below 0.2 M© yr -1 at 
£ = £ n = 6 Myr. 

Most sink particles move far away from the regions 
where they form; this is most likely a numerical effect aris¬ 
ing due to the timesteps and smoothing lengths of the sink 
particles being too large for proper orbit integration within 
the clusters. However, one can still identify five or six clus¬ 
ters (Figure [H middle panel), distributed asymmetrically 
within the CMZ. The most massive such cluster is located 
at {A, Y} ~ {50,100} pc (visual estimate; Figure 01 right 
panel). By measuring the cumulative mass within a given 
distance from the centre of the cluster, I estimate its size 
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Figure 8. Projected gas column density of models Base2 (left panel), Base3 (second panel), Base4 (third panel) and Base5 (last panel) 
at t = 6 Myr. The gas morphologies differ significantly, with no well-defined rings in any of the models. 
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Figure 9. Fragmentation rates in the five Base models. The al¬ 
gorithm calculates the change in total sink particle mass, so the 
actual rates may be slightly higher if many sink particles are ac¬ 
creted by the central SMBH. All models except Base3 show qual¬ 
itatively similar behaviour - initially small fragmentation rates 
lead to a surge once the CMZ ring begins fragmenting in earnest, 
at 3.5 — 4.5 Myr. In the Base3 model, fragmentation rates stay 
roughly the same throughout the simulation. 

to be R c \ ~ 10 pc and mass 2.8 x 10 6 M 0 . A more rigor¬ 
ous estimate of cluster size is not reasonable due to lack of 
numerical resolution and appropriate star formation physics 
within the model. Other clusters have sizes ~ 5 pc or less 
and masses ~ 0.7 — 1.6 x 10 6 M 0 . The typical star for¬ 
mation efficiency leads to star cluster mass estimates of 
1.4 x 10 4 — 2.8 x 10 5 M 0 , very similar to the masses of 
the observed clusters. The orbits of the sink particle clus¬ 
ters, much like the gas orbits, are non-circular and may be 
strongly irregular. This i s consistent with observations of the 
motion of both A rches (jStolte et al.ll2008h and Quintuplet 
(jStolte et al.ll2014h clusters. 

The gas and sink particles are distributed almost sym¬ 
metrically around the Z = 0 plane. In Figure 0 I show the 
gas column density in a side view of the CMZ. Small varia¬ 
tions in gas height and position are visible across the CMZ, 
but the displacement from Z — 0 does not exceed the ring 
thickness h < 40 pc. 


The other four Base models exhibit very different gas 
morphologies in the XY plane. Figure [8] shows the gas distri¬ 
bution in those models at t = 6 Myr. All four gas distribu¬ 
tions are clumpy, but the large-scale gas distribution differs. 
Model Base2 (left panel) has a full ring, model Base5 (right 
panel) shows approximately two thirds of a ring with a large 
gap in the bottom-left direction, model Base4 has prominent 
spiral streamers which may coalesce into a ring structure in 
several Myr. Finally, model Base3 has a very weak (low den¬ 
sity contrast) and highly off-centre elliptical structure, with 
very dense gas at the apocentre in the positive-X direction. 
It is interesting to note that even though the morphologies 
of the five models differ significantly, the formal ratio of 
semimajor to semiminor axis lengths, when defined as the 
maximum ratio of mean diameters of the dense gas features 
in two perpendicular directions, is rather similar, ranging 
from 1.34 to 1.53 (Table [fl last column). 

The morphological differences are reflected in the frag¬ 
mentation rates (Figure [9] and Table [I] columns 4-6). In 
models Base4 and Base5, several prominent bursts of frag¬ 
mentation occur earlier than in Basel, starting at f ^ 
3.4 — 4.0 Myr. The peak fragmentation rates in these mod¬ 
els are slightly higher than those in Basel, while the mean 
rates over the significant fragmentation period are similar. 
In model Base2, significant fragmentation starts at a similar 
time as in Basel, but proceeds slower. In model Base3 the 
fragmentation rate stays approximately constant, between 
5 — 10 M 0 yr -1 , throughout the simulation. The total mass 
of sink particles varies from 2.7 x 1O 7 M 0 in Basel and Base2 
to 3.8 x 1O 7 M 0 in Base4 and Base5, with an average value 
(3.2db0.5) x 1O 7 M 0 . The vertical structure of all four models 
is very similar to that of the Basel. 


4-2.3 Mixing of gas and angular momentum 

A general outcome of the five Base simulations is that an ini¬ 
tially smooth and stable distribution of gas becomes clumpy 
and anisotropic. This can be understood as a consequence of 
mixing of gas and its angular momentum. During the AGN 
phase, a large amount of gas (almost the whole CMZ) is 
pushed together into a narrow ring. This both increases the 
surface density of this gas, making it prone to gravitational 
instabilities, and efficiently mixes the angular momentum of 
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Figure 10. Distribution of CMZ gas angular momentum in the 
model Basel at t = 0, 1.5, 3 and 6 Myr (black solid, red dashed, 
blue dot-dashed and green long-dashed lines, respectively). The 
gas initially has a flat angular momentum distribution, which 
becomes progressively more squashed toward a single peak at 
1.5 X 10 4 pc km/s. 

gas, ensuring that the ring disperses slowly after the AGN 
switches off. 

To show this, I plot the angular momentum distribution 
of gas at t = 0, t = 1 Myr (end of AGN phase), t — 3 Myr 
and t — 6 Myr in Figure [TO] Initially, the angular momentum 
distribution is flat throughout the disc, since each annulus 
contains the same mass of material. AGN feedback pushes 
gas with low angular momentum (central parts of the disc) 
to merge with the high angular momentum gas, leading to 
a bump at middle values l ~ 10 4 pc km/s; the outer edge of 
the disc is not perturbed radially at t = 1 Myr. Later, this 
process continues, and by t = 3 Myr, the bump has moved 
to larger values due to mixing with the material at the out¬ 
skirts of the disc. Even by the end of the simulation, ma¬ 
terial is still being pushed together by spiral perturbations 
and gravitational instabilities; there is, however, a growing 
high angular momentum tail, which arises due to viscous 
spreading of the ring material. The total mass of material 
with very low angular momentum increases after the AGN 
switches off; this represents clumps of gas falling inward and 
coming very close to the SMBH, potentially feeding it (see 
Section 15.41) . 

4.3 Shorter AGN episode models: t q = 0.7 Myr 

The five Base models consistently produce gas rings, or parts 
of rings, that are too large and form sink particle clusters too 
late as compared with observations of the CMZ structure. 
Motivated by this, I run four models with t q = 0.7 Myr, 
labelled Short 1 to Short4 (see Table [I] for details and main 
results). As before, the models have stochastically different 
initial conditions, leading to different fragmentation histo¬ 
ries and gas morphologies at the end of the simulation. 

Qualitatively, the evolution of the Short models is very 
similar to that of the Base models. Initially, the AGN wind 
pushes disc gas together into a narrow ring, which becomes 
gravitationally unstable and fragments over several dynam- 



Radius [pc] 
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Figure 12. Radial density profiles of CMZ gas in model Shortl. 
Top panel: t = 0 (black solid), t = 0.5 Myr (red dashed) and 
t — 1 Myr (blue dot-dashed). Grey points are a sample of SPH 
particle densities at t = 1 Myr. Bottom panel: t = 1 Myr (black 
solid), t = 3 Myr (red dashed) and t = 6 Myr (blue dot-dashed). 
Grey points are a sample of SPH particle densities at t = 6 Myr. 

ical times. The shorter AGN episode duration results in a 
smaller initial extent of the ring, which is mirrored by the fi¬ 
nal sizes as well (see Figures fill fl2l and fl3l) . The ring eventu¬ 
ally has a radial extent from ~ 50 to ~ 80 pc from the centre, 
comparable to the observed CMZ size. The density contrast 
on each edge is not as sharp as in the Base models, but still 
noticeable, so that the ring can be identified as such (Figure 
1121 blue dashed line and points in the bottom panel). The 
ring evolves faster: since the expansion of the ring happens 
at constant velocity, and the subsequent evolution depends 
on the dynamical time, which increases linearly with radius, 

One WOuld expect tfrag,short /^frag ,base — ^q,sh°rt/^q,b ase — 0.7. 

In fact, rapid fragmentation begins 2.5 — 3 Myr after the 
start of the simulation (Figure fl4]) . almost 30% earlier than 
in the Base models, validating this simple explanation. The 
star clusters forming during this period would have ages 
consistent with observations of the Arches and Quintuplet 
clusters. Significant fragmentation has finished in all mod¬ 
els by 6 Myr, so the expected star formation rate today is 
< 0.2M© yr -1 . The masses and sizes of the star clusters are 
more difficult to determine than in the Base models, since 
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Figure 11. Projected gas column density of models Short 1 (left panel), Short2 (second panel), Short3 (third panel) and Short4 (last 
panel) at t = 6 Myr. The perturbed gas rings are smaller than in the Base models, with a variety of structures present. 



Direction [deg] 

Figure 13. Average radius (black solid line), half-diameter (red 
dashed line) and displacement (blue dot-dashed line) of the dense 
gas ring, defined as n > 100 cm -3 , as function of direction coun¬ 
terclockwise from left in model Short4 at t = 6 Myr. The radius 
varies between ~ 70 and ~ 110 pc. The half-diameter varies be¬ 
tween ~ 75 and ~ 105 pc, and the SMBH is displaced from the 
centre of the ring by about 10 pc. 


the older sink particle clusters have lost more members due 
to numerical effects. The total masses of sink particles form¬ 
ing in the four simulations are 2.4—2.9x 10 7 Mq, only slightly 
lower than in the Base models. 

The final shape of the gas rings varies among the four 
Short models (Figure fill), but the rings are generally more 
regular than in the Base simulations. Their ellipticity is sim¬ 
ilar to the longer-duration runs, a/b — 1.27 — 1.55 (Tabled 
last column), but the displacement of the SMBH from the 
centre generally does not exceed 10 pc. As an example, I plot 
the radius, half-diameter and displacement of the gas ring in 
the Short4 model in Figure fl3l The ring half-diameter varies 
between 75 and 105 pc, quite similar to the 60 and 100 pc 
of the model bv iMolinari et all (2011). Sgr A*is displaced 
from th e centre of t he ring by ~ 10 pc, less than in the 
IMolinari et al.l ( 2011 ) model, but that value is less strongly 
constrained than the size of the ring itself. 

Overall, the Short models give a final gas distribution 
in better agreement with the observed shape and size of the 



Time [Myr] 

Figure 14. Gas fragmentation rate as function of time in the 
four Short models. The bursts of significant fragmentation occur 
at t = 2.5 — 3 Myr, earlier than in the Base models. 

CMZ than the Base models. However, the stochastic evolu¬ 
tion of the gas ring results in large uncertainties regarding 
various more detailed parameters (such as the number, size 
and mass of star forming regions and young clusters), pre¬ 
venting a more rigorous constraint on past AGN activity 
from being made. I now turn to the four models designed to 
show the effect of varying the AGN duration and luminosity 
more significantly. 


4.4 Parameter variation 

Figure fl5l presents the gas and sink particle morphology of 
the final four models with varying parameters. From left to 
right, these models are: L0.3, where the AGN luminosity 
is decreased to 0.3LEdd, TO.3, where t q — 0.3 Myr, T1.5, 
where tq — 1.5 Myr, and L0.667T1.5, which has both L = 
0.667Z/Edd and tq — 1.5 Myr. 

Both L0.3 and TO.3 models show similar morphology. 
The time evolution of the models is also similar - a short or 
weak burst of AGN activity only affects a small inner part of 
the CMZ before switching off. The disturbance propagates 
outwards, but grows progressively weaker and is unable to 
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Figure 15. Gas column density of models L0.3 (left panel), TO.3 (second panel), T1.5 (third panel) and L0.667T1.5 (last panel) at 
t = 6 Myr. The gas morphologies differ significantly, with various structures visible. Note the different scale in the third panel. 
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Figure 16. Fragmentation rate in models L0.3 (black solid line), 
TO.3 (red dashed), T1.5 (blue dot-dashed) and L0.667T1.5 (green 
long-dashed). In all models, fragmentation is stronger soon after, 
or even during, the AGN phase, gradually decreasing with time 
as the gas distribution relaxes. 

produce gravitational instability in the outer parts of the 
CMZ. Some clumps form in the central regions; they produce 
sink particles and create spiral waves, which affect the whole 
disc and are clearly visible in the morphology plots. The TO.3 
model produces more sink particles (total mass at t = 6 Myr 
is -Mto.3 — 2.3 x 10 7 Mq and A/lo.3 — 1.0 x 10 7 M©), 
but in both models, essentially all fragmentation occurs at 
t < 2 Myr (see Figure fl6l black solid and red dashed lines). 
It is interesting that model TO.3 has a single burst of frag¬ 
mentation at £ ~ 1.6 Myr, which occurs when the perturbed 
material falls back to the centre of the disk and reaches 
highest density. 

In model T1.5, the prolonged action of the AGN un¬ 
binds most of the CMZ gas and allows it to disperse over 
a very wide region (Figure [15] third panel; note that the 
physical scale is different than in the other two). This can 
be expected, because even though most of the outflow en¬ 
ergy escapes in directions perpendicular to the CMZ, some 
of that energy is retained and adds to the pure momentum 
push in unbinding the gas disc. The gas density in the cen¬ 
tral regions is reduced and there is very little gravitational 


fragmentation (Figure fl6l blue dot-dashed line). Some spiral 
perturbations are present, and a ring of relatively dense gas 
forms at a distance d ~ 100 pc from the centre, but this ring 
contains only ~ 5 x 10 6 M© of gas. Sink particles disperse 
throughout the volume, with roughly one third of the total 
A/'t 1.5 = 1.3 X 10 7 M© escaping from the central 200 pc of 
the Galaxy. 

The model L0.667T1.5 provides the same energy input 
to the CMZ as the Base models, but spread more gradually. 
Its evolution is markedly different from the Base models. A 
more gradual energy input results in a much lower density 
contrast in the CMZ, which then leads to more gas being 
expelled from the system and far less fragmentation (Fig¬ 
ure □m green triple-dot-dashed line). The total sink particle 
mass at t = 6 Myr is Mlo.667ti.5 = 1.3 x 10 7 M©. This 
mass is concentrated in just a few sink particles, which are 
distributed with no particular structure throughout the sim¬ 
ulation volume (Figure fl5l right panel). 

All four models with parameter variation exhibit very 
little fragmentation after t — 2 Myr, making the clusters 
formed in these simulations older than those observed in the 
Galactic Centre. In addition, the morphology of dense gas 
differs significantly from the observed ring. These discrepan¬ 
cies suggest that the set of parameters chosen for the Base 
and Short models - a 0.7 — 1 Myr long Eddington-limited 
burst of AGN activity 6 Myr ago - are the best fit to the 
Sgr A* activity that could have produced the dense CMZ 
ring and massive star clusters. 


4.5 Summary of results 

A brief summary of the main results of the simulations is as 
follows: 

• A burst of Eddington-limited AGN feedback lasting for 
0.7—1 Myr is able to transform an initially smooth gas disc 
in the Galactic centre into a clumpy asymmetric distribu¬ 
tion. 

• The morphology of the distribution depends strongly 
on stochastic processes - spiral and gravitational instabili¬ 
ties. It is not guaranteed, but seemingly quite likely, that 
an elliptical ring of dense gas forms in the CMZ, and that 
the major-to-minor axis ratio of the ring, as well as the ring 
size, is very similar to the observed one. 
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• Fragmentation of the CMZ into sink particles takes off 
several dynamical times after the AGN activity finishes. The 
dynamical time, and hence the age of the sink particle clus¬ 
ters, depends on the duration of the AGN episode. For the 
trmq = 1 Myr episode, fragmentation occurs ~ 3.5 —5.5 Myr 
after the start of the simulation, i.e. 0.5 — 2.5 Myr ago. If 
the AGN episode lasts for only t q — 0.7 Myr, fragmentation 
starts earlier, 3.5 Myr ago. The average fragmentation rate 
during the last 6 Myr is < 10M© yr -1 , with peaks reach¬ 
ing > 30M© yr -1 . Given a typical star formation efficiency 
on the mass scale of the sink particles, e* = 2 — 10%, the 
star formation rate in the CMZ is M* < 0.2 — 1M© yr -1 , 
with bursts up to ~ 3M© yr -1 . By t = 6 Myr, fragmen¬ 
tation in the Short models has dropped to < 2M© yr -1 , 
corresponding to a star formation rate of < 0.2M© yr -1 ; 
however, fragmentation is still ongoing in some Base mod¬ 
els. 

• The mass of sink particles formed in the Base simula¬ 
tions is (2.7 — 3.8) x 10 r M©, predominantly distributed in a 
few clumps, each with mass in the range (0.5 — 3) x 10 6 M©. 
Using the same star formation efficiency as above, the star 
clusters formed should have masses in the range 10 4 — 3 x 
10 5 M©. The Short simulations form a slightly smaller to- 


The rate of formation of new sink particles at the end of 
the simulation, when multiplied by the typical star forma¬ 
tion efficiency, produces a star formation rate estimate simi¬ 
lar to the observed value ~ 0.06 —0.08M© yr -1 . This value is 
at first sight puzzling on two accounts. It is hi gher than cou ld 
be ex pected from steady-state arguments ( Zubovas et al.l 


>20131 ). i.e. the current SFR would exhaust the molecular gas 
reservoir in much less than the age of the Galaxy. On the 
other hand, it is several orders of magnitu de lower than ex¬ 
pecte d from the surface density of the gas (|Longmore et al.l 
120131 5 . Explanations for the first mismatch include bar- 
induced inflow of material, wh ich replenishes the gas con - 
tent in the Galactic centre (|Rodriguez-Fernandez fc Combesl 
l2008i ). while the second can be explained by an unusually 
high level of turbulence in the Galactic centre molecular 
clouds, which r educes t he star format ion rate at a given 
surface density (|Rathborne et al.l 120141 5 . The AGN-induced 
fragmentation model can explain both phenomena simulta¬ 
neously. The elevated star formation rate is a transient phe¬ 
nomenon, caused by the AGN episode compressing the gas, 
and will likely drop by at least a factor of a few in the upcom¬ 
ing several Myr. Meanwhile, the molecular clouds forming 
from the fragmenting CMZ disc are subject to very high ex- 


tal mass of sink particles, but the cluster masses are more 
difficult to determine due to numerical effects. 

ternal pressure from the expanding AGN outflow (the Fermi 
bubbles: see Zubovas et al. 2013|j. which cause them to have 

• An AGN outburst that either has a lower luminos¬ 

higher densitv and turbulence (Elmeg 

men & Efremov 

1997), 

ity, L = 0.3LEdd, or a significantly shorter duration, t q = 

leading to a higher SFR threshold ( 

Longmore et al. 

20li; 


0.3 Myr, produces a perturbation to the CMZ that is too 
weak to create a dense ring with embedded young star clus¬ 
ters. Conversely, a longer AGN burst duration, t q = 1.5 Myr, 
disrupts the CMZ, leaving only a minor fraction of the gas 
bound in the Galactic centre. 

• An AGN outburst with the same total energy output 
as in the Base models, but spread on a longer timescale 
produces a weaker density contrast in the CMZ disc, which 
leads to an almost non-existent dense gas ring and far less 
sink particle formation. This type of outburst might still 
enhance the existing density contrasts in the CMZ, but is 
unable to shape it to the present state from smooth initial 
conditions. 


5 DISCUSSION 

5.1 Current properties of the CMZ 

The numerical models presented here adequately explain the 
origin of some of the peculiar CMZ properties. Notably, an 
Eddington-limited AGN episode which started 6 Myr ago 
and lasted for 0.7-1 Myr could have shaped an initially axi¬ 
ally symmetric CMZ disc into an elliptical, off-centre dense 
gas ring with embedded massive molecular clouds and star 
clusters. The total mass of dense gas and sink particles 
(which represent both stars and gas denser than resolvable 
in the simulation) is very similar to the observed value, al¬ 
though this is potentially an imprint of the initial conditions, 
and a more (less) massive initial gas disc would produce a 
more (less) massive dense gas ring at the end. The fact that 
such complex structures can be created from axisymmetric 
initial conditions suggests that lopsided nuclear structures 
should be a common feature of galaxies (see Section 15.31 be¬ 
low). 


The AGN-induced fragmentation model naturally ex¬ 
plains why the massive clusters in the CMZ are younger 
than the nuclear star cluster: the NSC formed from the gas 
that fed the AGN outburst, while the CMZ clusters formed 
from the perturbed CMZ gas. The ages of the sink particle 
clusters forming in the Short models are ~ 2.5 — 3.5 Myr 
(Figure [9]), consistent with observations of the Arches and 
Quintuplet clusters in the Galactic Centre. The precise ages 
of star clusters are impossible to determine from the sim¬ 
ulations, both due to insufficient resolution to resolve star 
formation, and due to the stochastic evolution of the gas 
ring, which leads to a variation in sink particle cluster ages 
among the simulations with identical parameters. 

Recent observations reveal several massive stars in the 
CMZ which do not belo ng to either of the young mas¬ 
sive clusters (|Pong et al.l 2014: 1. These stars have veloci¬ 
ties relative to their surrounding gas ranging from ~ 10 to 
~ 80 km/s; some of them might have formed in-situ, while 
others might have been ejected from parent clusters. In my 
models, both processes occur - some sink particles form in 
the dense gas filaments, representing clusters or associations 
too small to be resolved in the simulations; other sink par¬ 
ticles form in clusters and then quickly leave them due to 
three-body interactions and tidal stripping. In order to bet¬ 
ter understand the relative importance of these processes, 
I plot, in Figure 1171 the histogram of relative velocity be¬ 
tween sink particles and surrounding gas (defined as the 40 
nearest SPH particles, with contributions weighted by the 
SPH kernel) for the model Basel at t = t n = 6 Myr. A clear 
peak at u re i = 50 km/s is visible, with ~ 26% of the relative 
velocities in the range {13; 78} km/s (vertical dashed lines 
in the plot), corresponding to t he m inimum and maximum 
relative velocities observed bv IPong et al.l ( 20141 ). There is 
a large tail of high-relative-velocity sink particles, indicat- 
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Figure 17. Histogram of relative velocity between sink particles 
and surrounding gas in model Basel. A strong peak around v re \ = 
45 km/s is present, with a power-law tail of particles escaping 
from newly-formed clusters. The vertical lines mar k the range o f 
relative velocities, 13 to 78 km/s, observed bv lDong et al.l d2014fl . 


ing that a lot of sink particles are ejected from clusters, 
as one would not expect stars formed in-situ to have rela¬ 
tive velocities larger than the background velocity dispersion 
g — 150 km/s. This may be partly a numerical effect due 
to insufficient detail in the star formation process in the 
model, as some sink particles might in reality be gas blobs 
which escape the clusters, disperse and do not form stars. 
Nevertheless, I predict that stars with significantly higher 
velocities relative to the surrounding gas should exist in the 
Galactic Centre than have been discovered so far. 

Certain properties of the CMZ are not reproduced 
by this model. One such feature is the twist of the ring 
(Mol inari et al .1120 111 ). which is not seen in any of the mod¬ 
els presented here. This suggests that the twist is either an 
imprint of different initial conditions, created by a differ¬ 
ent process than AGN-induced compression and fragmenta¬ 
tion, or a much rarer effect of AGN activity than the effects 
that have been seen. A barred potential leads to vertical as 
well as radial oscillations of orbits and so could perhaps cre- 
ate the twisted shape required by the first two explanations 
([Binnev et al.|[l99lf : lHasan et alll 19931 ). 

Another feature which the models presented here do not 
capture is the dynamical properties of the star clusters, i.e. 
their sizes, density distributions and mass functions. The 
lack of these properties is due to the crude star formation 
prescription, which only follows the initial collapse to a self- 
gravitating high-density core, and due to insufficient mass 
resolution in the models. We may expect a high integrated 
star formation efficiency in the cores due to the high external 
pressure preventing the parent molecular clouds from dis - 
persing (|Elmegreen &; Efremovl Il997l : IZubovas et al.l 120141 ) 
and a slightly top - heavy IMF due to the aforementioned 
high pressure (|McKee fe Tanll2002h and the gas dynamics 
being similar t o the collect-and-collap se model of induced 
star formation (| Whitworth et al .11 19941 ). 
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5.2 Footprint of past Sgr A* activity 

Comparison of the results of the Base, Short and variation 
models allows me to put constraints on the past activity of 
Sgr A*. The AGN episode 6 Myr ago cannot have been much 
weaker than Eddington-limited and cannot have lasted much 
less, or much longer, than 0.7-1 Myr. These results are con¬ 
sistent with the cons traints pro vided by the exist e nce of the 
Ferm i bubbles ([Zubovas et al.l 1 20 111 : IZubovas &; Navakshinl 
2012) and the presen ce of a 6 Myr old nuclear star cluster 
([Paumard et alJl2006h , which is eviden ce of a star formation 
event and a likely AGN feeding event (jHobbs &; Navakshinl 
l2009h . Together, these pieces of evidence strongly support 
the argument that an Eddington-limited AGN episode oc- 
cured in the Galactic centre ~ 6 Myr ago and lasted for 
~ 0.7 — 1 Myr. The dynamical footprints of this evidence 
can still be seen today, providing a glimpse into the past 
state of the Galaxy. Similar evidence in other galaxies may 
help investigate their recent evolution (see Section 15. 3D . 

ISofuel ( 20031 ) suggested that star formation in the Galac¬ 
tic centre might proceed in recurrent episodes, with foot¬ 
prints of many such episodes visible throughout the Galaxy. 
In his model, the North Polar Spur was formed from a star- 
burst that occurred 15 Myr ago, and the Galactic Centre 
Lobe is the remnant of a several-Myr old starburst. My re¬ 
sults are consistent with this picture. A burst of AGN ac¬ 
tivity shuts off the feeding reservoir for a period of several 
Myr, but allows for more feeding later (see Section 15.41 be¬ 
low); this process may cause AGN luminosity to fluctuate 
in short bright bursts separated by periods of insignificant 
activity. Since each AGN episode may lead to compression 
of surrounding gas and a period of enhanced star formation, 
the star formation history in the Galactic centre should be 
composed of short bursts separa ted b y longer periods of in¬ 
activity, exactly as suggested bv ISofuel ( 20031 ). 


5.3 General properties nuclear rings affected by 
AGN feedback 


Many galactic centres contain l arge amounts of d ense 
gas a rranged in ri n gs or discs (iKormendv &; Kennicuttl 
2004; iKnaoenl 120051 : iKrips et al.l 120071 : iBoker et al.l 20081 ). 


driven inflow of material from kpc scales ( Combes & Gerin 

1985 

: lYuan & Kud 1997; Sakamoto et al. '19991: Combes 

2003 

; Kormendy & Kennicuttl 

2004 Sheth et al.l 120051: 

Rodriguez-Fernandez & Combes 

20081), although there are 


|2004 : ISirchenko &; Moiseevl 12006 ). A bar-driven gas in¬ 
flow process predicts several morphological features of the 
rings : ellipticity (since the rings follow triaxial X 2 or- 

bits:|Morris &; Serabvnll 19961 : iKnapen et al.ll200l : iKim et all 

l20nl ). more massive gas clumps close to t he o rbit cross¬ 

ing (since material moves slowest there; iBenedict et al.l 
1 19961 ) and an azimuthal gradient of star cluster ages (since 
star formation is tr i ggere d at the gas inje ction points; 
IFalcon-Barroso et al.l 120071 : IBoker et al.l 120081 ) . All three 
properties have be en o bserved in som e galaxies, but are 
not ubiquitous (|Knapenl 120051 : IFalcon-Barroso et al.1120071 : 
08 |). 


IBoker et al.l [200. 

The model of AGN-induced compression of the CMZ 
gas presented in this paper might appear as an alternative 
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to the bar-driven inflow model. However, it is more appro¬ 
priate to view the two models as complementary. The trans¬ 
port of gas into the central parts of the galaxy is mediated 
by bar potential or caused by minor mergers. Afterward, the 
nuclear gas ring is affected by both the aspherical gravita¬ 
tional potential and episodic AGN feedback. The latter has 
several important effects on such nuclear gas rings: facilitat¬ 
ing the formation of more chaotic ring morphology, causing 
bursty and elevated star formation, and producing a corre¬ 
lation between star cluster ages and ring sizes. 

5.3.1 Ring morphology 

An AGN outflow disturbs the gas ring radially and in¬ 
duces instabilities. These alter the shape of the ring, mak¬ 
ing it irregular, rather than axially symmetric. This effect 
has been observed both directly, in NGC 613, where AGN 
activity is associated with disturbances in gas morphology 
([Falcon-Barroso et al.ll2014h . and can be inferred for many 
other gal axies, where irregular morphologies o f gas rings are 
common ([Sarzi et al.1120071 : iBoker et al.ll2008h . 

The location of most massive molecular clouds is 
also random in the AGN-induced compression model. 
In our Galactic Centre, the star forming regions Sgr 
B and Sgr C are located clos e to the ends of the 
CMZ ellipse ([Molinari et alJfjoTl ), as predicted by a bar- 
driven inflow model. Howe ver, other ga l axies do not show 
such regular structures ( Knapenl 1 20051 : I Sarzi et ahl 120071 : 
iFalcon-Barroso et alJ 12014 ). so a more stochastic process of 
generating molecular clouds is preferred there. 

Finally, AGN-induced fragmentation also leads to ran¬ 
dom locations of star formation and hence to random az¬ 
imuthal distribution of star cluster ages. This is called a 
“popcorn model”, as opposed to the “ pearls on a string” 
model predicted by bar-driven inflows (iBoker et al.l l2QQ8h . 
In this latter model, star and cluster formation occurs pre¬ 
dominantly at the ends of the ellipsoidal ring, where the 
gas from the galactic bar enters the ring; the young clusters 
then travel along the ring, producing an azimuthal gradient 
of stellar ages, with the youngest stars at the ends of the 
ellipse and oldest stars just upstream of those points. Many 
galaxies show azimuthal age gradients of star clu sters i n nu¬ 
clear rings, but this isn ’t a ubiquitous feature (| Sarzi et al.l 
120071: iBoker et, al.l 120081) . 

5.3.2 Star formation rates 

The mass flow rate of a bar-driven gas inflow varie s little ove r 
time, since the stellar b ar is typically long-lived ([Kim et al.l 
l201ll : ISeo &; K im 2013). This results in an approximately 
constant star formation rate and predicts a rather uniform 
spread of star cluster ages. Observations, however, reveal 
that star formation in rings is mo re likely to be bursty 
([Allard et al.ll2006l : ISarzi et ahll2007l j. suggesting a need for 
episodic perturbations. AGN feedback can provide just such 
perturbations, increasing the SFR over the long-term aver¬ 
age. AGN feedback also explains why the observed star for¬ 
mation efficiency in nuclear sta rbursts is much higher than 
in typical star-forming regions ([Kennicu ttl ll998h : AGN out- 
flows create conditions of very high pressure around the gas 
ring, whi ch prevent gas from escaping and enhance star fo r- 
mation (Elme gree n & Efrem ovl ll997l : Zub ovas et al. 120141 "). 


5.3.3 Ages of star clusters 

The sequence of events at the basis of my model - AGN 
fuelling and outburst, mixing of gas in the nuclear ring and 
formation of molecular clouds and star clusters - predicts a 
clear difference between the age of star clusters in the nuclear 
ring and the time since the last significant SMBH accretion 
episode, which should coincide with the last star formation 
episode in the central few parsecs of the galaxy. Therefore, 
the age difference between the youngest stellar population in 
the nuclear star cluster (NSC) and the youngest star clusters 
in the nuclear gas ring can be expressed as 


At a g e — ATdyn (Rring) CK -Rring/<^gal5 


( 7 ) 


where R r ing is the radius of the nuclear gas ring and cr ga i is 
the velocity dispersion in the galaxy spheroid. K is a con¬ 
stant, here assumed to be independent of ring properties 
or the larger galactic environment; this assumption is sup¬ 
ported by the results of the Base and Short models (see 
Section Oil . Using the parameters of the Milky Way, I find 
K rsj 3. Current observations cannot easily distinguish the 
youngest populations in NSCs, since those populations prob¬ 
ably comprise only a small fraction of the total NSC mass 
(in our Galactic centre, the mass of stars formed in the burst 
6 Myr ago is a few per cent of the total NSC mass), but in 
the future, it should be possible to test this correlation. 

One example where this correlation can be seen is 
th e galaxy NGC 16 IT w hich hosts a nuclear star burst 
([Alonso-Herrero et al.ll200ll : iKonig et al.ll2013l j. The galaxy 
contains a very young NSC, t ag e,NSC = (7 ± 1) x lO 6 yr 
(|de Griis et al l2005h : its central v elocity dispersion is cr ga i 
164 ±8 km/s (Hinz & Rieke 200 6|). The nuclear gas ring ex¬ 


tends radially from 100 to 350 pc |Xu et al.ll2014t h with most 
star-forming clumps found at a radius iAing = 230 pc. Us¬ 
ing equation © above and K — 3 gives At age ,NGCi6i4 = 
4.2 Myr, i.e. the star clusters in the nuclear gas ring should 
be ~ 2.8 Myr old. The precise ages of the star clusters in the 
nuclear ring are not known, but the presence of large HII re¬ 
gions and evidence of ongoing star form ation suggests tha t 
they are younger than 5 Myr ([Alonso-Herrero et al.lE002 b 
consistent with the model prediction. In fact, there is evi¬ 
dence that the n uclear starburst in NGC 1614 is propagating 
outward ([Alonso-Her rero et al.ll200lh . which may simply be 
a result of longer dynamical time, and hence slower evolution 
of a perturbed gas ring, at larger radii. This explanation pre¬ 
dicts an apparent outward propagation of the starburst with 
velocity v ou t — &gai/K ~ 55 km/s, in excellent agr eement 
with observations (^ ou t,obs = 60 km/s [Alonso- Her rero et ahl 
2(1(11 ). 

The fact that a burst of nuclear star formation follows 
the AGN activity episode with a delay gives an opportunity 
to put constraints on possible past AGN activity in a galaxy. 
This was discussed for our Galaxy in Section 15.21 above. 
An interesting example is M82, the prototypical starburst 
galaxy. It contains a large number of young star clusters i n 
its ~ 250 pc core , with ages ~ 10 Myr ([Satvapal et al.|[l997l : 
iMelo et ahl l2005h . The velo city dispersion in t he bulge of 
M82 is a ~ 100 =b 20 km/s ([Gaffney et al.|[l993h . Assuming 
that the core is a ring-like structure, the time since the most 
recent AGN burst should be 17.5 Myr or so - long enough 
for the AGN to have switched off and its most direct effects 
to have faded, but short enough for an outflow to remain 
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visible llKing et, al.lboi l'). A part of the outflow observed in 
M82 may be a fossil AGN outflow, rather than purely due 
to starburst activity. On the other hand, if one assumes that 
the star formati on prop a gates inwards from the larger-scale 
older starburst (|de Griisl]200ll h the time delay between the 
young star cluster ages and AGN outburst implies a powerful 
AGN activity episode starting 5 d= 2.5 Myr ago; one might 
expect obvious evidence of recent powerful AGN activity 
in this case, which is not det ected (o nly a low -luminosity 
AGN is detected in M82; iMatsumoto fe Tsurulll99^ b. The 
outward propagation of the starburst is also supported by 
observations of young s tar cluster ages, giving a propaga¬ 
tion velocity ~ 50 km/s (|Satvapal et al.lfl99 7), qualitatively 
consistent with the model of an AGN-induced starburst. 


5.4 AGN fuelling 

Although this is not the primary goal of the paper, the re¬ 
sults presented herein offer several suggestions regarding fu¬ 
elling of AGN in the presence of feedback. 

The usual mechanism invoked for transporting gas from 
outer regions of galaxies toward the centre is gravitational 
torques from stellar bars (ICombesI 120031 . also see Section 
15.31 above). However, they tend to cause gas to accumu¬ 
late in nuclear rings several hundred pc in size. Further 
gas infall is t hwarted by the bar re s onances and change 
of or bit types (jRegan k, Teubenll2003l : iGarcfa-Burillo et al.l 
2005). One possible mechanism for bringing gas closer to 
the AGN is by action o f nested bars (IShlosman et al.l Il989: 
iMacieiewski et ahll2002h . but only a few galaxies show evi¬ 
dence for such structures. Other explanations invoke tran¬ 
sient effects, such as increased gas turbulence due to star- 
burst activity in the central kiloparsec, which enhances mix¬ 
ing of gas of different angular momentum and allows gas 
to pe netrate further into the galactic centre dHobbs et al.l 

imj. 

My results show that mixing of gas angular momentum 
can be achieved by the AGN itself, by means of feedback 
acting upon gas of spatially varying angular momentum (see 
Figure [TO] and Section [4. 2. 31) . In the case of the CMZ, this 
mixing is not very efficient, because gas is already in ordered 
motion around the centre, i.e. the mean angular momentum 
of the gas is large and the formal circularization radius of the 
whole CMZ is also large. This allows only a small fraction of 
all the gas to lose most of its angular momentum via the ac¬ 
tion of local gravitational torques produced by gravitational 
instabilities in the CMZ. In a more general case, where gas 
is falling in through the bulge, the motion is generally less 
ordered and the total angular momentum of the gas is small. 
In that case, mixing of gas, enhanced by AGN feedback, can 
cause gas clouds with very different angular momenta to 
collide and fall toward the centre. This s cenar io was inves¬ 
tigated analytically by iDehnen fe Kind (12013 ), who found 
that SMBHs can easily be fed at their Eddington limit by 
repeated short bursts of feedback. The results of my mod¬ 
els, showing efficient angular momentum mixing, are consis¬ 
tent with the basic picture of this work. In particular, this 
irregular mixing might lead to massive reservoirs of gas ac¬ 
cumulating very close to the SMBH, within the sphere of 
influence. This might be the case in Andromeda, where the 
double nucleus shows a massive reservoir of gas within a few 


pc of the SMBH, but w ith little star formation (|Tremainel 
Il99fil: iBender et, al.ll200. lr J) . 

5.5 Outstanding issues 

The study presented here is a first step in understanding the 
effects of AGN activity upon nuclear gas rings. As a result, a 
number of simplifying assumptions had to be made in order 
to isolate the salient effects. Below, I discuss the importance 
of the effects neglected by making these assumptions. 

5.5. 1 Gas equation of state 

The models presented in this paper adopt a realistic gas 
heating-cooling prescription above T = 10 4 K, but prevent 
gas from cooling below this value. This is done because cool¬ 
ing above 10 4 K is less dependent upon detailed effects of 
gas self-shielding, which my models cannot resolve. In princi¬ 
ple, the temperature floor stabilizes the disc against gravita¬ 
tional instabilities and washes out possible small-scale struc¬ 
tures in the perturbed gas rings, leading to a lower fragmen¬ 
tation rate than might be expected if the gas was allowed 
to cool to lower temperatures. Since the density threshold 
for fragmentation varies as nj oc T 3 (eq. [5]), a temperature 
floor at T = 100 K would reduce this threshold by a fac¬ 
tor 10 6 . However, the peaks in the gas density distribution 
are extremely sharp and tall, rising more than 6 orders of 
magnitude above the typical densities in the gas disc/ring. 
Therefore, even such a significant change in the temperature 
floor should not affect the fragmentation rate significantly. 

Feedback from star formation might also affect the frag¬ 
mentation rate, since stars can heat their surroundings and 
hamper or shut off further fragmentation in their vicinity 
(|Batel 120091 )). Proper implementation of feedback, however, 
requires much higher mass resolution than in the simulations 
presented here, sufficient to resolve massive stars and small 
stellar associations. 


5.5.2 Gravitational potential 

The set up of the models presented in this paper makes an 
assumption of a spherically symmetric back ground poten¬ 
tial. The presence of a Galactic bar (IBissantz fe Gerhardl 
120021 : iBabusiaux fe Gilmore! 120051 ) makes the actual poten¬ 
tial triaxial, which may significantly affect the results of 
these simulations. Qualitatively, an asymmetric potential 
should enhance the asymmetry in the final gas distribution, 
so the main findings of this paper should still be valid. 

5.5. 3 Initial conditions 

The initial gas distribution used in the simulations presented 
here is necessarily strongly idealised. One reason for choos¬ 
ing a smooth disc, rather than a realistic lumpy distribution, 
is that any enhancement of gas fragmentation is likely to be 
even stronger when more realistic gas distributions are con¬ 
sidered. This would occur because a realistic gas distribution 
has pre-existing clumps and density inhomogeneities, which 
would be enhanced by a passing AGN outflow. On the other 
hand, some parts of the disc may be cleared, rather than 
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compressed, so the induced star formation rate would be 
smaller in a more realistic situation. 

Given that most of the outflow impacting the CMZ is 
deflected perpendicularly, it might seem that the outflow 
has only a weak effect on the Galaxy on kpc scales. This 
was found in some previ ous simulations of high-r edshift, 
therefore gas-rich galaxies dGabor fe Bournaudll2014h . How¬ 
ever, one has to keep in mind that even though the out¬ 
flow is collimated, its opening angle is still la rge, of order 
7r/2 for each cone ([Zubovas &; Navakshinl 11201 2lj. This value 
is con sistent with recent observational results |Nardini et al.l 
l2015h . Therefore, a large part of the galaxy is affected by the 
outflow directly (radially), a nd the rest of the gas may be 
compressed perpend icularly (fsll k 2005; Gaibler et al.l [2012: 
IZubovas et al.ll2013h . If the initial conditions are made more 
realistic, the outflow might escape in even more directions, 
leading to an even larger effect on the host galaxy. 


5.5.4 Numerical resolution 

Another important aspect of the simulation is the numerical 
resolution. With the resolution chosen for the simulations, 
the smoothing length of CMZ gas particles, /isph, is between 
0.1—2 times the scale height of the disc at the beggining of 
the simulation, with values /isph > Hcmz reached only on 
the surface of the disc. This means that the vertical thickness 
of the CMZ disc is typically resolved with a few tens of 
particles, making the vertical structure marginally resolved, 
and hence fr agmentation in the disc should not be artificially 
suppressed dNclsoif2000 b 

The masses of fragments forming in the CMZ depend 
only weakly on numerical resolution. While initially they 
form from single SPH particles, they quickly accrete and/or 
merge with other particles, growing in mass to several hun¬ 
dred time the SPH particle mass. Therefore, the masses of 
those fragments are generally resolved, from a purely nu¬ 
merical perspective. Their masses would undoubtedly be 
affected if additional physics (such as stellar feedback) is 
added to the simulation. 

Numerical resolution also affects the interaction be¬ 
tween the AGN wind and the CMZ gas. The momentum and 
energy carried by each virtual particle is scaled to the total 
SPH particle number, so this direct interaction should not 
be strongly affected. However, the effective boundary of the 
CMZ disc, where gas can absorb the wind energy and evap¬ 
orate vertically, is always ~ Hcmz thick, since virtual parti¬ 
cles are typically absorbed in a region of this length scale. 
Therefore, with higher resolution, this boundary becomes 
physically thinner and less dense. This leads to less evapora¬ 
tion from the disc, but more importantly, the gas is heated 
to higher temperatures by the interaction with the AGN 
wind. Therefore, it cools more slowly (since t coo \ oc T 1 / 2 
for Bremsstrahlung cooling, dominant at the large temper¬ 
atures relevant in this case) and both evaporates and passes 
energy to the surrounding matter more rapidly. As a result, 
the disc might expand more both in the radial and vertical 
directions. Given that this effect depends not just on res¬ 
olution, but also on the gas heating-cooling prescription, a 
thorough analysis is beyond the scope of this paper. 


6 CONCLUSION 

I presented the results of hydrodynamical simulations of a 
brief AGN activity episode and the associated outflow affect¬ 
ing the gas disc around the Galactic centre. The simulation 
set up is designed to mimic the properties of our Galaxy, 
and the time since the AGN outburst, 6 Myr, is fixed by 
the age of the young stars in the central parsec. The initial 
distribution of the gas has the size and mass of the Central 
Molecular Zone, but is distributed in a smooth disc. 

With this setup, I show that an Eddington-limited AGN 
outburst lasting for 0.7—1 Myr can by itself transform the 
smooth gas disc into a narrow, elliptical and off-centre ring 
of dense gas clumps and young star clusters. The masses of 
star-forming regions and star clusters agree rather well with 
the observed features in the CMZ - the Sgr B2 molecular 
cloud and the Arches and Quintuplet star clusters - and the 
shape of the ring is also similar to observations. However, the 
stochastic nature of gravitational instabilities, which cause 
the formation of structures, means that there is a lot of 
variation among the models and more concrete predictions 
are difficult to make. 

Several properties of the CMZ cannot be explained with 
this model. Most importantly, the CMZ is known to be 
twisted perpendicularly to its plane, a feature not repro¬ 
duced in any of my simulations. This twist may potentially 
be explained by the action of an asymmetric background po¬ 
tential, i.e. the Milky Way bar, but further work is required 
to test this hypothesis. Certain other observed features of 
the CMZ are consistent with a purely bar-driven gas in¬ 
flow model. A recent AGN outburst and associated outflow 
should be considered as an addition to this model, not a 
replacement. In the future, I plan to model the response 
of a more realistic barred galactic environment to an AGN 
outflow. 

Given that past AGN activity can explain many fea¬ 
tures of the CMZ, these features can be used to constrain 
the duration and luminosity of the AGN episode. These con¬ 
straints agree with previous constraints based on the Ferm i 
bubbles ([Zubovas et al.l l201ll : IZubovas &; Navakshinl l2Q12h . 
Together, these aspects are dynamical footprints of past 
AGN activity, and point to the possibility of investigating 
activity history in other galaxies as well, by finding appro¬ 
priate dynamical structure that might have been affected by 
past AGN outflows. 

Nuclear gas rings and star formation within them is a 
common feature in many galaxies. While bar-driven gas in¬ 
flows can explain many features of these rings, some aspects, 
such as high star formation efficiencies during repeated star- 
bursts, disturbed morphologies and radial gradients of star 
cluster ages, can be more easily explained by the action of 
AGN outbursts. Overall, AGN outbursts have a significant 
influence upon the properties of host galaxy centres, and 
future observations will help quantify this influence. 
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